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We present measurements of the magnetic susceptibility, heat capacity and electrical re- 
sistivity of Pui_ r Lu^Pd3, with x=0, 0.1, 0.2, 0.5, 0.8 and 1. PuPd3 is an antiferromagnetic 
heavy fermion compound with = 24 K. With increasing Lu doping, both the Kondo and 
RKKY interaction strengths fall, as judged by the Sommerfeld coefficient y and Neel tem- 
perature Tff. Fits to a crystal field model of the resistivity also support these conclusions. The 
paramagnetic effective moment /j e ff increases with Lu dilution, indicating a decrease in the 
Kondo screening. In the highly dilute limit, /j e ff approaches the value predicted by interme- 
diate coupling calculations. In conjunction with an observed Schottky peak at ~60 K in the 
magnetic heat capacity, corresponding to a crystal field splitting of ~12 meV, a mean-field 
intermediate coupling model with nearest neighbour interactions has been developed. 

I. INTRODUCTION antiferromagnetic, with a transition temperature 

~24 K, and a G-type structure, where near- 
The AnPd 3 series of compounds, with est neighbour moments are aligned antiparallel. 
An=U, Np, or Pu, are rare examples of actinide ^ same stud y found that the hi § h temperature 
intermetallic compounds in which the 5/ elec- resistivity shows a Kondo-like behaviour, in- 
terns are well localised around the ionic sites, creasing with decreasing temperature. This be- 
UPd 3 is a very interesting compound which ex- haviour, together with a high Sommerfeld coef- 
hibits four quadrupolar ordered phases below ficient deduced from recent heat capacity mea- 
8Ki, whilst there are indications that NpPd 3 surements 5 , led us to make a further investiga- 
may also show quadrupolar order at low tern- ^ on or " me properties of PuPd 3 . 
peratures^. These two compounds crystallise ur aim has been to study how the compe- 
in the double-hexagonal close-packed (dhcp) tition between the Ko ndo effect and the RKKY 
structure, in contrast to PuPd 3 which adopts exchange interaction affects the physical prop- 
the AuCu 3 structure, with lattice parameter erties of this compound . Th i s may be accom- 
a =4.105A. This reflects the increasing locali- plished by doping with a non -magnetic ion, 
sation of the 5/ electrons as shown by recent which i ncrea ses the distance between localised 
photoelectron spectroscopy measurements^. /-electrons, and hence decreases the RKKY in- 
Early measurements of the bulk properties teraction. Finally, the single-ion properties of 
and neutron diffraction studies^ show it to be the Pu ion may be investigated in the highly di- 
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FIG. 1 : (Color online) Lattice parameters and 
transition temperatures for Pui T Lu x Pd3. 



lute limit. 

In this work, we present in section HI] the ex- 
perimental details and in section [HI] the mea- 
surements of the magnetic susceptibility, and 
heat capacity of Pui x Lu x Pd3. These mea- 
surements are then analysed using a localised 
moment mean field model in section [IVj Fi- 
nally, section [V] presents measurements of the 
electrical resistivity and Hall coefficient of 
Pui_ A Lu Y Pd3, which are assessed in terms of a 
simple crystal field model. 



II. EXPERIMENTAL DETAILS 

Poly crystalline samples of PuPd3, LuPd3 and 
Pui^LUx-Pds, with x = 0.1, 0.2, 0.5, and 0.8, 
were produced at ITU by arc melting appropri- 
ate amounts of the constituent elements under a 
high purity argon atmosphere on a water-cooled 
copper hearth, using a Zr getter. The AuCu3 
structure was confirmed by x-ray diffraction 
for each sample, and the lattice parameters are 
shown in figure \T\ The data show a linear de- 
pendence of the lattice parameter with increas- 
ing Lu dilution, in accordance with Vegard's 
Law, which also confirms the stoichiometry of 
the samples. In addition, there also appears to 
be a linear dependence of the transition tem- 
perature, Tjy, with doping. These temperatures 



were determined from magnetic susceptibility 
and heat capacity measurements described in 
the next section. At there is a maximum 
in the susceptibility, %, and hence this temper- 
ature was determined by numerically differen- 
tiating the data to find -k = 0. In the heat ca- 
pacity, C p , there is a lambda step at 7V, which 
was determined by differentiating the data to 

find the minima of The values of 75v de- 
duced from these two measurements are in close 
agreement, whereas the inflexion points of the 
resistivity data do not correlate with TV as de- 
termined from % or C p . Nevertheless, the resis- 
tivity inflexion points show the same decreasing 
trend with Lu-doping as 7^. Errors in Tn quoted 
in table U were determined by the width in tem- 
perature of the lambda step for C p or the step in 

d% 
dT- 



The magnetisation and susceptibility were 
measured using a SQUID magnetometer (Quan- 
tum Design MPMS-7), whilst the heat capac- 
ity was determined by the hybrid adiabatic re- 
laxation method in a Quantum Design PPMS- 
9 for PuPd 3 and LuPd 3 , and in a PPMS-14 
for Pui_ x Lu x Pd3. Small samples with mass 
less than 5 mg were used for the heat capacity 
measurements so that the decay heat does not 
significantly affect the measurements. The X- 
ray diffraction, magnetisation and heat capacity 
measurements were made immediately after the 
preparation of samples in order to minimise the 
effects of radiation damage. 



Finally, thin parallel- sided samples of each 
composition were extracted, polished and 
mounted for electrical transport measurements. 
As these measurements were made some three 
months after synthesis, we observed significant 
radiation damage which manifested in a high 
residual resistivity at low temperatures. This 
prompted us to anneal the samples at 800° C for 
12 h, and to remeasure the electrical transport 
properties. The remeasured data is presented in 
section El 
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FIG. 2: (Color online) The magnetisation at 
2 K of Pui x Lu x Pd3. Lines are guide to the eye. 



III. MAGNETISATION AND HEAT 
CAPACITY MEASUREMENTS 

A. Magnetisation 

Figure [2] shows the magnetisation at 2 K, 
which shows that the Pu-rich compositions are 
not saturated at 7 T. This is not surprising be- 
cause we expect a J = 5/2 ion to be saturated at 
a field > 100 T, when the splitting between the 
lowest two CF levels is ^>2 K. The magnetisa- 
tion of Puo.2Luo.8Pd3 shows some evidence of 
saturation, however. The magnetic susceptibil- 
ity is shown in figure [31 and the inverse suscep- 
tibility in figure HI The data in the paramagnetic 
phase above the Neel temperature are well fitted 
by a modified Curie- Weiss Law 
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where N is the number of Pu atoms in the com- 
pound, and Gcw is the paramagnetic Curie tem- 
perature. We recall that in the Weiss mean field 
theory, (-Gcw) Ocw corresponds to the (anti-) 
ferromagnetic transition temperature. However, 
this theory does not take into account single ion 
effects such as the crystal field which are ex- 
pected to be significant in PuPd3. 

There is a field dependent residual sus- 
ceptibility %o(H) which mainly arises from 



impurities, the encapsulation and the sam- 
ple holder. In addition, the Pauli sus- 
ceptibility of the conduction electrons may 
also contribute to %o an d can be estimated 
from the electronic Sommerfeld coefficient 
of LuPd 3 , 7+^ =3.2(1) mJmol -1 K- 2 , which 

yields x Pau ii « 5.2(2) x 10~Vs T 1 f.u. -1 . This 
is significantly lower than the observed values 
of the residual susceptibility, which are of the 
order of 10 3 jUgT 1 f.u. , indicating that the 
conduction electron susceptibility contribution 
is negligible. The fitted parameters to the Curie- 
Weiss relation for each sample are given in ta- 
ble HI The quoted error is deduced from calcu- 
lating the covariance matrix of the parameters 
from the covariance matrix of the data 6 assum- 
ing that this is diagonal and proportional to £ 2 
where £ is the standard error in the measured 
moment or heat capacity as determined by the 
MultiVu software supplied by Quantum Design. 
Figure @] shows the inverse susceptibility of the 
different compositions with the residual suscep- 
tibility %o subtracted. 

The magnitudes of the effective moments 
are all significantly higher than the LS-coupling 
value, 0.85 Any crystal field interaction will 
only decrease this effective moment because as 
the crystal field split levels become further sepa- 
rated and hence thermally de-occupied, their an- 
gular momentum will cease to contribute to the 
moment. The effective moment with zero crys- 
tal field splitting in intermediate coupling on the 
other hand is approximately 1.4 /jb, as calcu- 
lated using the theory outlined in section [IVl 
This value may be decreased slightly by a large 
crystal field, and suggests that we should use in- 
termediate coupling to calculate the single-ion 
properties of Pu 3+ . 

An alternative reason for the higher than ex- 
pected effective moment may be due to some 
high moment paramagnetic impurity. How- 
ever, analysis of the x-ray diffraction patterns 
showed that the only detectable impurity is 
LU2O3 which is non-magnetic. There may also 
be some trace amounts of oxides of Pu which is 
not observed in the diffraction pattern. PuC>2 is 
a Van Vleck paramagnet 7 which may contribute 



3 





. PUPClg 




PU 0.9 LU 0.1 Pd 3" 




\ A PU 0.8 LU C.2 Pd 3- 
\ □ PU 0.5 LU 0.5 Pd 3 
V - PU 0.2 LU 0.8 P( V 




s^-a^V + LuPd 3 


^^^^ 

-^^+ ++ . , , , 


^ s * S * aaa ^^A g □ „ 

' '•* , °5oaoaoSo 


" +++ ++++ + + + + + + 





10 20 30 40 50 

T(K) 

FIG. 3: (Color online) The magnetic 
susceptibility of Pui x Lu x Pd3. 
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FIG. 4: (Color online) The inverse magnetic 
susceptibility of Pui x Lu x Pd3. Thin solid lines 
are fits to the modified Curie-Weiss Law for 
each compound. The mean field intermediate 
coupling calculation is shown as a thick solid 
line, whilst the dashed line shows the single ion 
susceptibility calculated from the crystal field 
in LS-coupling. 



to the impurity term in equationQ], whilst PU2O3 
is an antiferromagnet with an effective moment 
of 2.1 figr. However, one would need approx- 
imately 6 mol % PU2O3 in PuPd3 order for it 
be responsible for the increased effective mo- 
ment compared to the LS-coupling expectation, 
at which concentration it should be detectable 
in the X-ray diffraction pattern, which is not the 



case. Furthermore, the enhanced value of /u e ff 
for PU2O3 which also has a free ion Pu 3+ con- 
figuration suggests that intermediate coupling is 
appropriate in these cases. 

B. Heat Capacity 

The heat capacity at zero field is shown in 
figure |51 whilst details of the results in applied 
fields up to 14 T are in figure [61 We note that 
at high temperatures, C p tends to the classi- 
cal Dulong-Petit limit, 3NR = 99.8 Jmol l K 1 . 
For Puo.9Luo.1Pd3, the derivative in the heat ca- 
pacity shows two minima, which stem from the 
step-like nature of the transition. The higher 
temperature inflexion point at ~22 K corre- 
sponds well with the peak in the inverse suscep- 
tibility, but the lower temperature peak at ~21 K 
does not match any feature in the magnetic sus- 
ceptibility. Nevertheless, these two anomalies 
raise the possibility that there may indeed be 
two transitions in this compound. Moreover, 
as can be seen in Figure [61 the heat capacity 
of Puo.8Luo.2Pd3 also shows indications of two 
transitions. 

An estimate of the electronic specific heat 
C e i = yT and Debye temperature dp was ob- 
tained using the approximation 

c ^ + _^y ( 2) 

which is valid at low temperatures (T <C Go), 
from a plot of C p /T vs T 2 shown as an inset in 
figure \5\ However, the magnetic heat capacity 
complicates the determination of y because the 
Neel temperature is very low in some of the Lu- 
rich compounds. This increases the low tem- 
perature C p and hence the estimate of y from 
the straight line intercept. For this reason, for 
the Lu-rich compositions, we show the results 
of fitting the data in the region above 7Jy (y + ) in 
addition to that below 8 K (y~) in table H For 
Pu-rich compositions, the data above will 
be affected by the Schottky anomaly at approx- 
imately 17 K, and will be unreliable. Thus it 
appears from these estimates that the electronic 



4 



Composition 


T N (K) 


Y^mJmor'K- 2 ) 


6d (K) 

T>T N 


ftff to Pu l ) 


Ocw (K) 


c P 


X 


T>T N 


T <8K 


PuPd 3 


24.4(3) 


24.9(4) 




75(1) 


306(19) 


1.02(1) 


-39.9(1) 


Puo.9Luo.1Pd3 


22.7(4) 


22.5(1) 




133(8) 


304(18) 


0.97(2) 
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1.15(2) 
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7(2) 




99(49) 
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1.25(1) 
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TABLE I: Transition temperatures and other parameters derived from the magnetic susceptibility 

and heat capacity of Pui x Lu T Pd3. 
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FIG. 5: (Color online) Heat capacity of 
Pui A Lu x Pd3 at zero applied magnetic field. 
The inset shows C p /T vs T 2 and solid lines 
there show fits to C/T ~ y+AT 2 . LuPd 3 data 
is shown in orange in the inset 



heat capacity first increases with increasing x 
until x ps 0.5, whereupon it falls as x rises fur- 
ther. The spread in the fitted parameters when 
data from different ranges of temperatures in the 
region 25 < T < 40 K for y+, and 2 < T < 8 K 
for y~ was taken as an estimate of the errors in 
these parameters. 

The Debye temperature 0£> was determined 
from fitting the high temperature data, and ap- 
pears to be independent of Lu-doping. This sug- 
gests that the phonon contribution to the heat ca- 
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FIG. 6: (Color online) Heat capacity of 
Pui Y Lu x Pd3 in applied magnetic field. The 
arrows indicate the direction of increasing field. 



pacity is constant through the series. A good es- 
timate of this contribution is given by the heat 
capacity of the non-magnetic isostructural com- 
pound LuPd3, which also has a negligible elec- 



5 











PU 0.9 LU 0.1 Pd 3 




PU 0.8 LU 0.2 Pd 3 




PU 0.5 LU 0.5 Pd 3 




" PU 0.2 LU 0.8 Pd 3 




" — CF Calculation 











Temperature (K) 

FIG. 7: (Color online) Deduced magnetic and 
electronic contribution to the heat capacity 

Cp ag = C f (P Ul _ x Lu x Pd 3 ) - (C p (LuPd 3 )-y+7:). 
The solid line is a crystal field calculation. 
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FIG. 8: (Color online) Deduced magnetic 
contribution to the entropy of Pui x Lu x Pd 3 , 
calculated by numerically integrating C™ g /T. 



tronic heat capacity, Ylu = 3.5(5) mJmol K . 
We have thus extracted the additional electronic 
and magnetic heat capacity of Pui x Lu x Pd 3 by 
subtracting that of LuPd3, as 

C^s = C p (P Ul _ x Lu x Pd 3 ) - (C p (LuPd 3 ) - y+ T) 

(3) 

This extracted quantity, scaled by the Pu con- 
centration, is shown in figure [7J 

The magnetic heat capacity for all the com- 
pounds shows a peak at ~60 K, which we at- 
tribute to a Schottky anomaly from the crys- 
tal field (CF) splitting. The cubic CF on 
the Pu 3+ ions splits the six-fold ground mul- 
tiplet (J = | in LS-coupling) into a doublet 
and quartet. The energy gap, A CF , between 
these two levels determines the temperature of 
the Schottky peak, such that A CF ~12 meV 
corresponds to a peak at 60 K. The magni- 
tude of this peak, however, is determined by 
whether the doublet (C p =6.3 Jmol _1 K _1 ) or 
quartet (C p =2 Jmol *K *) is the ground state. 
The data in figure [7] thus suggest a doublet 
ground state. 

This is supported by the magnetic entropy, 
shown in figure [8l deduced by numerically in- 
tegrating the magnetic heat capacity, S(T) = 



T C mag 

f 2 -y-dT . From the very low heat capacity of 
PuPd 3 at low temperatures, we believe the mag- 
netic entropy from to 2 K is negligble, and 
have not included this range in the integration. 
The value of the entropy at the Neel tempera- 
ture is approximately R In (3) for PuPd 3 , which 
is above the value, R\n{2), expected for a dou- 
blet ground state. If the electronic heat capac- 
ity Yp u = 76 mJmol -1 K -2 is subtracted from 
the integral, then we obtain S(Tn) ~ i?ln(2.4). 
The remaining discrepancy may be due to (i) an 
incomplete subtraction of the phonon contribu- 
tion, as the heat capacity of LuPd 3 may not be 
exactly analogous to the phonon heat capacity 
of PuPd 3 , and (ii) a larger value of y (see the 
discussion in section UVT). 

Finally, the inset to figure [5] shows what ap- 
pears to be a hump at ~ 17 K in the heat capacity 
of PuPd 3 . This feature was initially attributed to 
the Schottky anomaly from the CF splitting in 
reference 5 . However it is more likely due to a 
Schottky anomaly from the splitting of the dou- 
blet ground state in the ordered phase, and in- 
deed such a feature is observed in the mean field 
calculations in the next section. A splitting of 
A MF = 3.5 meV gives a peak at ~17 K, which 
is reasonable. 
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FIG. 9: (Color online) Calculated dependence 

of the CF parameter Lq and the nearest 
neighbour exchange parameter J7 on Lq subject 
to constraints described in the text. The vertical 
line indicates the optimal parameters. 



IV. MEAN FIELD CALCULATIONS 

As noted in section IIII A[ the measured ef- 
fective moment for both PuPd3 and the doped 
compounds is approximately 1 /Jg/Pu, in con- 
trast to the expected LS-coupling value of 
gj\/J(J+ 1) = 0.85 /jfi from a Hund's rule 
6 H 5/ / 2 ground state for a Pu 3+ ion. LS- 
coupling is a good approximation when both 
the Coulomb (#c) an d spin-orbit (^4 ) interac- 
tions are large but J-Cq ^> ^4o- hi contrast, when 
rtc <fC ^4o> JJ-coupling is a better approxima- 
tion, whereupon we obtain // e ff = 2.86 jig. In be- 
tween these limits, for the case of intermediate 
coupling, the effective moment is a function of 
#c an d rtsoi and the full Hamiltonian, including 
both these terms and the crystal field (#" c f) and 
Zeeman interactions (^z) must be calculated. 

The strength of Hq and 9{ so , parameterised 
by the Slater (F k ) and spin-orbit (^) integrals, 
is fixed by the atomic environment of the un- 
filled shell electrons. Thus in practice, inter- 
mediate coupling refers to the case where the 
value of F k and £, are determined either from 
ab-initio (Hartree-Fock) calculations, or from 
experimental measurements using optical Spec- 
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FIG. 10: (Color online) Measured and 
calculated magnetic heat capacity of PuPd3. 
The estimated electronic heat capacity, 

AC™ ag /L = p T p is given in the bottom 
panel and has a mean value of 
122 mJmol l K 2 . 



troscopy. Using parameters determined exper- 
imentally by Carnall from the spectra of dilute 
Pu 3+ in LaCls 9 , the effective moment is 1 .44 /ig. 
It is conceivable that in a metallic system like 
Pui_ x Lu T Pd3 there may be small changes to F k 
and ^ compared to the insulating salts on which 
the measurements o^ were made-™ Neverthe- 
less, a 10 % change in F k and ^ to make the sys- 
tem more LS-like only yields /j e ff =1.38 /jg. In 
order to obtain /j e ff ~1 /jg, we must double the 
magnitude of F k and £, compared to their mea- 
sured values, which is probably unphysical. 

Given that the crystal field interaction is 
small, as judged by the ~12 meV split between 
the doublet ground state and first excited quartet 
deduced from the heat capacity measurements, 
9{ c f has little effect on // e ff. Thus we believe 
that the lower than expected effective moment is 
most likely due to Kondo screening. Nonethe- 
less, a mean field calculation which can model 
the antiferromagnetic order and the single ion 
intermediate coupling behaviour is still valuable 
to interpret the heat capacity and magnetisation 
data. Such a calculation, carried out using the 
McPhase package- , is detailed below. 

We have assumed a nearest neighbour only 
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exchange interaction between 5/ electrons, 
which is reasonable given the G-type antiferro- 
magnetic structure where nearest neighbour mo- 
ments align in antiparallel. Thus, there are three 
free parameters in the calculations: two crystal 
field (CF) parameters, Lq and Lq, and one ex- 
change parameter J . There are two other non- 
zero CF parameters, but they are fixed by the 
cubic point symmetry of the Pu 3+ ions such that 
L\ = y/5j\4I% and h\ = -^TJlL%. It should 
be noted that the parameters used here corre- 
spond to the Wybourne normalisation 12 , rather 
than the usual Stevens normalisation (usually 
denoted B" 1 ). This is because the Stevens opera- 
tor equivalents Of are valid only within a single 
multiplet of given 7, whereas we now require 
operators that can span all the allowed J values. 

The two CF parameters are fixed by the re- 
quirement that they result in a doublet ground 
state with a quartet at ~12 meV. This fixes a 
relation between Lq and Lq as shown in fig- 
ure |9j The Neel temperature 7/v then fixes a 
relation between j? and the crystal field parame- 
ters, and finally the magnetisation below 7# was 
used to fixed all three values, yielding Lq = 320, 
L§ = 1250, and 3 = -0.204 meV. 

The magnetisation is calculated by including 
in the Hamiltonian a Zeeman term, — ^.g (L + 
2S) ■ H; numerically diagonalising the energy 
matrix and calculating the expectation value of 
the moment operator L + 2S. The calculated 
inverse susceptibility is shown as a solid black 
line in figure H] for comparison with the mea- 
sured data. Unfortunately, better agreement 
with the data within the constraints of the mean- 
field intermediate coupling model is only pos- 
sible by increasing the Coulomb or spin-orbit 
integrals to unphysical values. A more likely 
explanation is the suppression of the effective 
moment by Kondo screening, which is not con- 
sidered in the current model. 

The heat capacity is calculated by numeri- 
cally differentiating the internal energy, (U) — 
Y in E n exp(—E, l /kBT) /Z, by the temperature, 
and the entropy by subsequently numerically 
integrating this. The calculated heat capacity, 
shown in figure \\0[ shows a shoulder around 



~20 K in accordance with the data which arises 
from a Schottky peak due to the splitting of 
the ground state doublet in the ordered phase. 
We can also estimate the electronic heat capac- 
ity by subtracting this calculated from the 
measured electronic and magnetic heat capac- 
ity, C™ ag , the result of which is shown in the 
bottom panel of figure [TOl The spike near 7/v 
is due to the differences in the sharpness of 
the calculated and measured transitions in this 
temperature range. Overall however, the mean 
value, y =122 mJmol ! K 2 , over the full tem- 
perature range is in fair agreement with that de- 
rived from the low temperature heat capacity, 
y pu =76 mJmol _1 K -2 . 

Similar mean-field heat capacity calculations 
for the other compositions, where the exchange 
coupling J7 was reduced to reflect the lower T N , 
did not yield the broad transitions seen in fig- 
ure 13 but rather the sharp lambda anomalies ex- 
pected of an antiferromagnetic transition. Thus 
a subtraction to deduce their electronic specific 
heat becomes increasingly untenable. The broad 
transitions observed in the data are probably due 
to disorder in the system as a result of the Lu 
doping. 

Finally, the calculated internal fields in the 
model are 226 T (180 T) at 1 K (20 K), which 
agree well with the molecular field of 217 T 
determined from fitting the measured resistiv- 
ity using a simple CF model, as described in the 
next section. 

V. ELECTRICAL TRANSPORT 
MEASUREMENTS 

Electrical transport measurements were car- 
ried out using thin parallel-sided samples ex- 
tracted after crushing the poly crystalline buttons 
produced by arc-melting. The first transport 
measurements were completed some months af- 
ter the production of the samples, so there were 
significant aging effects in the Pu samples. This 
prompted us to anneal the PuPd3 sample, and 
re-measure its resistivity, resulting in a large de- 
crease in the residual resistivity po from 225 /u£l 
cm to 11 /u£l cm. Subsequently the resistivity 
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of the other compositions was also re-measured 
after annealing. The values of po show a rapid 
increase with Lu doping up to x = 0.5 there- 
after decreasing with x, as summarised in ta- 
ble HQ This is due to the increasing number of 
defects caused by Lu substitution. Finally, mea- 
surements of the Hall coefficient and longitu- 
dinal resistivity of PuPd3 and the x=0.1,0.2,0.5 
compositions in field were also carried out. 

The resistivity of LuPd3 is well fitted by the 
Bloch-Griineisen relation 



T r x x 

P = C-g- / - —dx (4) 

9^ Jo (e - 1)(1 —e x ) 

where x = T/Q D , with C = 181(1) jjQ. cm 
and D = 166(1) K. It was taken to be rep- 
resentative of the non-magnetic contribution 
to the resistivity of Pui Y Lu x Pd3, and used 
to estimate the magnetic resistivity as Ap = 
p(Pui -xhUxVdj)— p(LuPd3). This quantity is 
plotted in the case of zero applied magnetic field 
in figure [TU The in-field measurements showed 
little change from the zero field data, and the 
data for PuPd3 agree well with previous mea- 
surements 13 , albeit with a slightly lower resid- 
ual resistivity. 

Qualitatively, the behaviour of the resistivity 
may be divided into a high temperature Kondo- 
like regime, where the resistivity increases with 
decreasing temperature until ~50 K, followed 
by the onset of coherence, from where it falls 
sharply with temperature, and shows no clear 
anomaly at 7/y. At low temperatures, the re- 
sistivity follows an exponential temperature be- 
haviour, in contrast the power law behaviour ex- 
pected in metals from the Bloch-Griineisen re- 
lation. Electrons scattering from antiferromag- 
netic magnons will give rise to an exponential 
temperature dependence, as will spin-disorder 
scattering from the localised 5/ moments them- 
selves. A fit to the electron-magnon resistiv- 
ity^ assuming an isotropic magnon dispersion 
(0 = VA 2 +Dk 2 , such that^ 



p e _ ra = p + BA 2 ^e-^ T (5) 

3A 2 / A \ 2 " 
X + 2k B T + 15 \k B T J 

yields, however, a spin-gap A = 4 meV which 
is significantly larger than that expected from 
the heat capacity below Tpj, if the shoulder at 
17 K corresponds to a Schottky peak which 
arises from a gap of approximately 3.2 meV. In 
contrast, this very splitting between the doublet 
ground states in the ordered phase is predicted 
by the spin-disorder resistivity model described 
below. 

Above ~70 K, the resistivity is well fit- 
ted by a po — pilog(r) term 1 -, where po is 
the residual resistivity, and pi is proportional 
to the interaction between the conduction elec- 
trons and Kondo impurities. The fit is shown 
in figure QT1 with parameters po = 235.8(3) /u£l 
cm and pi = 22(1) /jQ. cm for PuPd3. The 
value of pi initially decreases with Lu doping 
to 15(2) fiQ, cm for x = 0.1 and 16(2) /jQ. cm for 
x = 0.2 but then increases to 27(3) /jQ. cm for 
x = 0.5. This increase suggests that the Kondo 
interaction is strengthened at half doping. 

The magnetic resistivity of Puo.2Luo.sPd3 
does not show the Kondo behaviour of the other 
compositions, but rather increases with increas- 
ing temperature with a plateau region around 
30-80 K. This behaviour and also the exponen- 
tial temperature dependence of the low temper- 
ature part of the resistivity is characteristic of 
a simple crystal field spin-disorder resistivity 
model 17 . This model is based on the scatter- 
ing of conduction electrons with spin s by a lo- 
calised moment J through an exchange interac- 
tion — 2Gs ■ J, giving the resistivity in the first 
Born approximation as 

3%Nm* ,. 7 l 

x E ( m ' S M\ s -J\ m sWi)Pifii' ( 6 ) 

m s ,m' s ,i,i' 



9 



300 



250 



200 



<3 150 

Q. 
< 

100 



50 



PuPd, 




FIG. 1 1 : (Color online) Zero-field magnetic 
resistivity of (Pu,Lu)Pd3 on a logarithmic scale. 
Inset is an enlarged view of the Puo.2Luo.8Pd3 
composition. Solid lines are fits to the relation 
po — pi log(r) described in the text. 



where the occupation factor for the crystal field 
level at E[ is p\ = exp( — ^) and the conduc- 
tion electron population factor is fy = 2[1 + 
— L )]~ . The wavefunctions and 



exp( 



k B T 



energies Ei are determined by diagonalising the 
crystal field Hamiltonian. 

In the absence of a crystal field, the 27 + 1 
degenerate spin orbit ground state levels yield a 
temperature independent resistivity given by 



3nNm* 
he 2 E F 



G 2 (g-l) 2 J(J + l) (7) 



In our case, with a J = | multiplet in a cu- 
bic crystal field, and in the absence of a mag- 
netic field, equation © reduces to a sum of 
exponential functions, because the wavefunc- 
tions |\|/;) are fixed, and the crystal field pa- 
rameter can only change the splitting A CF be- 
tween the quartet and doublet. There is thus a 
universal behaviour, with the resistivity tending 

to Ps_f at T 3> A CF , and then falling exponen- 
tially as the temperature falls below some level 
such that the excited crystal field states are no 
longer populated. This temperature is approxi- 
mately 0.4A CF , so the model suggests a splitting 



Composition 


Po 




B m f 


pic 
B mf 




(pD. cm) 


(pQ. cm) 


(T) 


(T) 


PuPd 3 


12.0(1) 


876(37) 


217(23) 


226 


Puo.9Luo.1Pd3 


75(1) 


118(5) 


88(4) 


149 


Puo.sLuo.2Pd3 


97(1) 


35(1) 


213(6) 


143 


Puo.5Luo.5Pd3 


172(1) 


53(2) 


92(5) 


129 


Puo.2Luo.8Pd3 


11.5(1) 


17.8(1) 


51(1) 


88 



TABLE II: Parameters for the spin-disorder 
resistivity model for Pui_ x Lu A Pd3 described in 
the text. The crystal field parameter 
54=0.041(1) meV was determined using the 
data for PuPd3, and thereafter fixed in the 
fitting of the other compositions. 



A CF « 3.2 meV, as the drop off in resistivity oc- 
curs around 15 K. This is significantly smaller 
than the splitting deduced from the Schottky 
peak at ~60 K but is similar to the splitting of 
the ground state doublet (A MF ) in the ordered 
phase, as determined by the shoulder at ~17 K 
in the heat capacity data. 

In order to accommodate this splitting, we 
introduce a molecular field B m f term. The low 
temperature exponential increase is then gov- 
erned primarily by the split doublet, with a sec- 
ond exponential step at higher temperatures due 
to the crystal field splitting. It turns out that 
this second step is not observed in the case of 
PuPd3 because the two steps merge into each 
other. Indeed a fit to the data below 7/v with all 
parameters in the spin-disorder resistivity model 
varying freely yields a CF splitting of 14 meV, 
in agreement with the heat capacity data. As 
B m f decreases in line with with increasing 
Lu doping, whilst the CF splitting remains con- 
stant, the two steps become more pronounced in 
the calculations. These two steps are observed 
in the case of Puo.2Luo.8Pd3, but the second step 
is masked by the Kondo screening in the other 
compositions. 

Using this simple crystal field model we ob- 
tained the parameters shown in table [TTJ Fig- 
ure [12] shows the resulting fit to the data with 
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FIG. 12: (Color online) Zero-field magnetic 

resistivity of (Pu,Lu)Pd3 with residual 
resistivity po subtracted. Solid lines are fits 
using the crystal field model described in the 
text 
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FIG. 13: (Color online) Temperature 
dependence of the Hall coefficient Rh = Pxy/B 
and magnetoresistance at 9 T (inset). The solid 
black line is a fit using the measured magnetic 

susceptibility as described in the text 



Po subtracted. The CF splitting was fixed for 
all Lu doped compositions to the value deter- 
mined from fitting the PuPd3 data. The molec- 
ular field determined from mean field calcu- 
lations, 5j£ f , with exchange parameters J = - 
0.202, -0.199, -0.1925, and -0.174 meV, for 
compositions x=0.1,0.2,0.5 and 0.8 respectively, 
is also shown in the table. Apart from the case 
of Puo.8Luo.2Pd3, the fitted is consistently 
lower than the calculated B 1 ^. This is due to the 
overestimation of TV in the mean field approxi- 
mation, because we have estimated J from TV, 
and is proportional to 3 . Thus, is also 
overestimated. 

Puo.8Luo.2Pd3 shows an upturn at low tem- 
peratures, which cannot be accounted for by 
the current model. In addition, this upturn af- 
fects the fit by decreasing the ratio between the 
maximum and minimum resistivity, and hence 

pf}f- The exponential increase in the resis- 
tivity also occurs at a higher temperature and 
over a broader temperature range in this com- 
position than in the others, which explains the 
anomalously high B m f. These features may be 
artefacts of the sample, because whilst an aged 
Puo.9Luo.1Pd3 showed the upturn at low tem- 
peratures, the annealed sample did not, whereas 



X = 





0.1 


0.2 


0.5 




32.9(4) 


12.1(2) 


1.5(1) 


0.9(1) 


Ri 


-6.97(8) 


-2.31(3) 


-0.74(1) 


-0.028(1) 



TABLE III: (Color online) Fitted Hall constants 
of Pui 1 L u \T J d3. All values are in 
(lO^m 3 ^ 1 ). 



both aged and annealed Puo.sLuo.2Pd3 samples 
showed the upturn. Furthermore the resistivity 
of the aged Puo.sLuo.2Pd3 sample is lower than 
that of the annealed sample. This suggests that 
the annealing had not fully repaired the radi- 
ation damage, and thus the resistivity may be 
strongly affected by crystallographic defects. 

Nevertheless when the fits were repeated us- 
ing the calculated as fixed parameters, the 

fitted p^ f changed by less than 10%. The fits 
thus showed that the /-conduction electron in- 
teraction decreases with Lu doping, with a slight 
increase for the x = 0.5 composition compared 
to x = 0.2 and x = 0.8, in agreement with the fits 
to the Kondo parameter p 1 . 

We now turn to the electrical transport prop- 
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erties in an applied magnetic field. The temper- 
ature dependence of the Hall coefficient Rh and 
the magnetoresistivity p xx in 9T are shown in 
figure [131 Pxx(T) follows the same behaviour as 
the zero field resistivity shown above, with the 
exception that there is a small peak just below 
the Neel temperature as shown in the inset to the 
figure. This was not observed previously and is 
reminiscent of the superzone scattering near Tm 
in the heavy rare earths^. 

The temperature dependence of the Hall ef- 
fect may be described phenomenologically by a 
scaling of the magnetisation, as in 

p xy (T)=R B+R m M(T) (8) 

where Rq is the ordinary and R\ the extraor- 
dinary Hall constant. In figure [131 the solid 
lines show fits of the Hall coefficient to this re- 
lation using the measured magnetisation data. 
Both the ordinary and extraordinary Hall con- 
stants were found to decrease with Lu doping, 
as shown in table [nil R\ is proportional to 
the conduction-/ electron exchange interaction 
strength G discussed above, so the decrease in 
its magnitude further indicates that this interac- 
tion becomes weaker with Lu doping. 

Finally, we observed that the magnetoresis- 
tance, p(B), is linear with applied field for all 
samples, and showed a negative slope at high 
temperatures and a positive slope at low tem- 
peratures, except for LuPd3 where the slope was 
always positive. We interpret the negative mag- 
netoresistance to be a sign of the Kondo effect at 
high temperatures, with the normal metallic be- 
haviour giving a positive magnetoresistance at 
low temperatures. 



VI. CONCLUSIONS 

We have completed extensive bulk properties 
measurements on antiferromagnetic PuPd3 and 
the pseudo-binary compounds Pui x Lu x Pd3. 
The transition temperature was found to de- 
crease linearly from 7]v =24.4(3) K in PuPd3 to 
7(2) Kin Puo.2Luo.sPd3. 



Heat capacity measurements show a Schot- 
tky anomaly at ~60 K, which was interpreted 
as arising from a crystal field splitting between a 
doublet ground state and an excited state quartet 
at ~12 meV. The deduced Sommerfeld coeffi- 
cient, y was found to be significantly higher than 
that expected for the free electron model, with a 
value of the order of 0. 1 Jmol K 2 determined 
by fitting the data directly and by subtracting the 
calculated magnetic and measured phonon con- 
tributions. Direct fits to the data suggest that y 
decreases with increasing Lu substitution. The 
magnetic heat capacity was calculated using a 
mean field model which showed that the shoul- 
der in the data corresponds to a splitting of the 
doublet ground state in the ordered phase with 
a gap of ~3.5 meV. The size of this gap is sup- 
ported by fits of the resistivity to a crystal field 
model. 

Magnetic susceptibility and magnetisation 
measurements also showed that the paramag- 
netic effective moment, jU e ff increases with Lu 
concentration, approaching the value expected 
in intermediate coupling. These observations 
arise from the Kondo interaction which sup- 
presses the effective magnetic moment but en- 
hances the electronic effective mass, m*. As 
the Kondo interaction decreases with Lu dop- 
ing, ^/ e ff is screened less, and yoc m* falls. 

Electrical transport measurements support 
this decrease in the Kondo interaction with in- 
creasing Lu concentration, x, as parameters pro- 
portional to the /-conduction electron coupling 
in a crystal field model of the resistivity and Hall 
effect were found to fall as x increases. 
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